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Reduction of total harmonic distortion in multilevel inverters is a difficult 
optimization problem that includes nonlinear transcendent equations having 
more than one local minima.The task becomes difficult as the level of the 
cascaded multilevel inverter increases.Optimization techniques helps in 
finding out the solution in a very less span of time giving the best possible 
results. This paper deals with the harmonic elimination of cascaded 
multilevel inverter with equal D.C. sources using a new optimization 
technique called Whale optimization technique. The objective of this paper is 
to find the best combination of switching angles to minimize the lower order 
harmonics and the total harmonic distortion is reduced. This aforesaid 
algorithm is applied to an 11-level cascaded H-bridge inverter.Results shows 
that Whale optimization algorithm gives better results and effectively 
performs the above mentioned task. 
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1. INTRODUCTION 

Multilevel inverters have gained a lot of interest in the past decade for increasing the voltage to 
different levels which make it suitable for medium and large power applications like statcom, electric drives, 
renewable energy, etc.The H bridge inverter previously produces an output voltage of Vdc, 0, - Vdc. This basic 
H bridge switching technique is extended to other circuits that can generate added output voltage levels. This 
multilevel output voltage gives a staircase waveform which is similar to sinusoidal waveform thus reducing 
the harmonic content in the output. There are different type of multilevel inverter topologies like neutral 
point clamped, flying capacitor and cascaded inverter out of which the cascaded multilevel inverter has 
gained a lot of popularity due to its improved quality and connection of independent d.c sources (SDCS) to 
each of the module so as to attain high power level at the output.[l]-[2]The cascaded multilevel 
inverter(CMLI) is very efficient in minimizing THD and gives better quality of power.lt is animportant 
topology as it is more simplein obtaining the preferred output voltage from several D.C. sources contrast to 
diode-clamped and flying capacitors type multilevel inverters. The utility of a multilevel inverter (MLI).[3] It 
has a modular structure with simple switching method and occupies a lesser space. By connectingadequate 
number of H-bridges in cascade along with a properly choosen switching scheme, a nearly sinusoidal output 
voltage waveform can be obtained. Each H Bridge operates with a different switching scheme which is used 
for harmonic control. If the number of H-bridge given by‘s’ then the no level of output voltage obtained per 
phase in CMLI is 2sh-1. Each H Bridge operates at a different delay angle resulting in a staircase waveform of 
the output phase voltage where the output voltage is the sum of all voltage generated by H Bridge.[4]-[5] The 
output phase and line voltage obtained and the efficiency of DC to AC conversion depends on the THD. 
Normally the output voltage waveform of a single phase inverter contains 33.5 percent of third harmonics, 20 
percent of fifth harmonics and 14.5 percent of seventh harmonics approximately. As the output voltage in 
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CMLI is obtained by the synthesis of many D.C. Sources to get the desirable a.c. value, hence it is essential 
to minimize the harmonics and to keep the total harmonic distortion (THD) within the specified limits 
prescribed by IEEE standards. [ 6 ]This paper is organized as follows. In section II the SHE and mathematical 
equations of SHE are given followed by the introduction of WOA is section III. In section IV WOA is 
applied in harmonic elimination problem to find out best switching angles is presented along with the results 
and discussion for anil-level inverter. This includes the calculation of THD, and the elimination of the 5th, 
7th, 11th, 13th harmonics at various modulation index is given followed by conclusion in section V. 




Eigure 2. Pulse width if CMLI 


2. SHE EQUATIONS FOR A CASCADE MULTILEVEL INVERTER 

In SHE-PWM technique the switching angles are generated to cancel a set of lower order 
harmonics at fundamental frequency[7]. Figure 1 shows configuration of CMLI. The output voltage 
waveform obtained in a CMLI is a staircase waveform as shown in Figure 2. As per the Fourier series of 
staircase waveform the even harmonics are less problematic with respect to the odd harmonics. Hence the 
harmonics those are to be eliminated are3, 5, 7..., up to k-1 harmonics. For the elimination of k-1 harmonics, 
k number of switching angles are to be generated.[8]-[9] In balanced three-phase, the third harmonics is 
negligible in line voltage whereas present in-phase voltage and therefore it is possible to reject 5*, 7th, 11th, 
etc from the line voltage waveform at a low switching frequency [10]. The output voltage waveform is given 
by 


^^(cos(na],) sinnwt) [ (1) 

(^n = l,3,5... "" J 

Where n is the harmonic number and constraints of 0<a^<a2 < a 2 <cc 4 <a^<^. The objective of this paper is 
to calculate the optimum switching angles 0 : 2 , 0 : 3 , 0 : 4 , 0:5 for an 11 level inverter having equal voltage 
sources i.e. Vi=V 2 =V 3 =V 4 =V 5 =V 

5 ?n = cos cti -I- cos 0:2 -I- cos a-^ -b cos 0:4 -I- cos 0:5 
0 = cos Sa-i -b cos 5a2 -b cos 50:3 -b cos 5 -b cos 5 

0 = cos la-^ -b cos la^ -b cos la-i, -b cos 7 -b cos 7 ( 2 ) 

0 = cos llui + cos lla2 -b cos 110:3 + cosll + cosll 
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0 = cos 13a^ + cos 13a2 -I- cos ISctg -I- cosl3 + cosl3 

Here Modulation index (Ma)=the ratio of the harmonic component of voltage to the maximum attainable 
voltage. If the maximum^' harmonic or fundamental voltage is Vi and the d.c voltage is equal to Vdc- [11] 

^ V^n _ cos a^+cos a2+cos 0:3+cos a:4+cos as 

^" 4 KVdc ~ k ^ ' 


/ V 2 

Again, the expression of THD= I( )where k=3, 5, 7, 11 ... (4) 

And Vi is the fundamental component of voltage and k is the harmonic number. 

2.1. Problem Formulation of Selective Harmonic Elimination Using WOA 

The objective function is given below [12]. For elimination of harmonics, each term of the 
function F(aj) should be equal to zero where Ojis the switching angles. 

F(ai) = (l]f=iCos(aj) - 5m)^ -I- (^^2f=iCos(5aj)) -I- 

(^Zf=iCos(7ai)) -h (:j^i;f=iCos(llai)) -h (5) 

(;j^2]f=]^cos(13a;)) With a constraints 0<aj< ^ 


3. WHALE OPTIMIZATION ALGORITHM (WOA) 

Whale optimization algorithm is anew metaheuristic optimization algorithm based on the hunting 
activities of humpback whales. [13]. The humpback whales are very brainy mammals that like to hunt krill or 
small fishes nearer to surface of the sea. They move in a 9 shaped path forming a unique bubble to deceive 
the small fishes.This method is quite unique in humpback whales. The mathematical model of WOA is 
divided into three parts encircling of prey followed by bubble net hunting method and Search the prey. 

3.1. Encircling Prey 

WOA technique considers thepresent best candidate solution is close to finest results. The 
leader among the humpback whale is defined first and other humpback whales followthe leader in 
hunting.The mathematical equation of above behaviour is given by 


D = |C.X * (t) -X(t)| 

(6) 

'X(t + l)=X *(t)-A-D 

(7) 

A — 2a.f — cT 

(8) 

r^i 

II 

NJ 

(9) 


Where‘t’ designates the present iteration X* is the location of the result attained which is updated 
after every reiteration for getting improved outcomes. Equation (6) gives the absolute value, ‘a’ shrinks 
linearly from 2 to 0 and ‘r’ is a random vector between 0 and 1. 

3.2. Bubble Net Attacking Method 

This method is the exploitation stage and the mathematical model of this method is as follows 

a. Shrinking Encircling Method 

This performance is attained by shrinking the value of ‘a’ in Equation (8).The function range of ‘A’ 
is by ‘a’. A will move alsodecreasedfrom (X, Y) to (X*Y*) for 0<A<1. 

b. Spiral updating position 

In hunting, humpback whales swim around the prey in shrinking circle and along a spiral-shaped 
path at the same time. Here 50% chance among shrinking, encircling or the spiral model to update the 
position of the whales and target to mimic the helix-shaped movement are as follows: 
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X(t-l-l) = D*e'’'cos27rl-l-A'*(t) (10) 

~D = |X‘(t)-^(t)l (11) 

and shows the distance of ith whale to the location of thetarget.b is a constant for describing the 
profile of logarithmic spiral, 1 is a random number [-1, 1] and is part by part multiplication. Humpback 
whales swim around the prey within a shrinking circle and along the spiral-shaped path simultaneously. 
Taking a chance of 50% for both mechanisms, we get equations as follows. 

X(t -I- 1) = D * e'’'cos27rl -f F • Dlf P>0.5 (12) 

X(t-l-l) =X *(t)-i'FifP<0.5 (13) 

Here P is a random number between [0, 1] 


3.3. Search for Prey 

To get the overall finest values updating has been done with randomly chosen search agent rather 
than the best agent. Xrand is the random whales in the current iteration. The equations are given as 

X(t + 1) (14) 

D=IC.F^-i| (15) 


X'.AXY X*.Y X.Y 



X* AX.Y*-AY X*. Y*-AY X Y* AY 

Figure 3. Bubble net search shrinking encircling 
mechanism 



Figure 4. Bubble net search spiral updating position 
mechanism 


3.4. WOA Algorithm and Flowchart [14] 

a State population size, parameters, maximum number of iterations 
b Set the repetition counter. 

c Fitness function F(aj) is considered and preliminary population is generated arbitrarily 
d Search agent are estimated and finest search agent is assigned the job 
e Repetition counter is increased t=tH-l 
f All constraints like a. A, C, P, etc are updated 

g Exploration and exploitation are applied as per the values of Pand A. 
h Repeat step 5, 6, 7 till end criteria satisfies, 
i Write the solution 
j End 
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4. RESULTS AND DISCUSSIONS 

Metaheuristic optimization algorithm plays an essential role in finding out optimum switching 
angles required for elimination of harmonics. As the level of the multilevel inverter increases, the no of 
variables in the objective function given in Equation 5 also increases which makes the problem more 
complicated. The whale optimization algorithm is a new algorithm which performs the above task 
effectively. This algorithm has many benefits as it can give all possible answer in the form of switching 
angles for any level of multilevel inverter at a quicker speed of convergence, etc. . In order to test the 
performance of algorithm, an eleven level inverter is taken and a supply voltage of 50V is applied to each of 
the Five H Bridge connected in cascade. In this paper, the maximum number of iterations taken for WOA is 
500, the population of search agents is 30 and the dimensions are 5. The program is run for many times and 
the best results were taken into consideration. Figure 5 shows the optimum switching angles generated at 
different modulation index with the help of above-mentioned algorithm and with the objective to minimize 
THD. It is observed in Figure 6 and Figure 7 thatTHD in line voltages decreases with the increase in the 
modulation index and the fundamental component increases with increase in modulation index respectively. 
The results show that the low-order harmonics of the line voltage are very lesser in magnitude and are 
suppressed effectively at various modulation index and results are given in table-1. Because of the 120 phase 
shifts between R, Yand B phase voltages, all third harmonics component in line voltages are also very small. 
The THD is found out to be 4.25% at a modulation index of 0.8. Higher order harmonics are of very less in 
magnitude in comparison to the lower order harmonics for various applications of electric drives, etc. It is 
observed that results are improved with WOA algorithm when it is compared with other algorithms. Figure 
8(a) and Figure 8(b) shows the line voltage and phase voltage waveform at optimum switching angles. At 
each modulation index, the program was run many times for finding out the optimum switching angles and 
the graph of the best score obtained with respect to modulation index is given in Figure 9. Finally the 
convergence curve in Fig-10 shows that the above algorithm converges faster and has the ability to escape 
local minima. A comparative analysis of the performance of whale optimization algorithm and other 
algorithms are also given in Table 2. 
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Table 1. THD Calculation of 5*^ 7th, 11th, 13th Harmonics 


Modulation 
index ’Ma’ 


Switching 

angles in 

radians 


- %THD 

a. 

a? 




0.5 

0.6201 

0.7935 

0.9957 

1.211 

1.4828 

1.04 

0.6 

0.4607 

0.76195 

0.9039 

1.0857 

1.2662 

1.1 

0.7 

0.1441 

0.4994 

0.7213 

0.9327 

1.2804 

0.7 

0.8 

0.116 

0.3295 

0.4764 

0.7887 

1.086 

0.6 



Figure 5. Switching Angles vs Modulation Index 



Figure 6. %THD vs modulation index 



Figure 7. Fundamental component vs modulation index 
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(a) 



(b) 

Figure 8. (a)Line voltage waveform (b)Phase Voltage Waveform 



Figure 9. Best score vs modulation index 



Figure 10. Convergence curve 
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Table 2, Comparision of THD Using Different Algorithm 


Optimization 

Modulation 


Switching angles 

in degrees 


■ %THD 

techniques[15][16] 

index ’Ma’ 

a. 

a? 




Newton -Raphson 
method 

0.8 

6.56 

18.93 

27.17 

45.12 

61.87 

5.55 

LSHADE /Epsin 
Technique 

not given 

5.49 

16.68 

28.58 

42.06 

59.46 

7.26 

WOA 

0.8 

6.6 

18.8 

27.3 

45.1 

62.2 

4.25 

Genetic Algorithm 

0.65 

3.54 

10.2 

14.66 

24.9 

38.3 

4.39 

Differential Evolution 
Algorithm 

0.7 

17.84 

39.1 

54.30 

63.62 

86.13 

5.23 


5. CONCLUSION 

In this paper,WOA algorithm is successfully applied to the problem of switching angle optimization 
in a multilevel inverter fed from equal sources.The proposed WOA algorithm scrutinizes for all probable set 
of results to contribute the least THD. From the results, we draw the conclusion that the suggested approach 
is able to reduce the lower order harmonics while retaining the fundamental output voltage.The5*, 7th, 11th, 
13this found out to be less than less than 1% whereas the higher order odd harmonics up to39th was found 
out to be 4.25%. It is observed that with WOA algorithm the results are comparatively better when compared 
with other algorithms along with its ability to escape local minima. 
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